Thermally induced coherence in a Mott insulator of bosonic atoms 
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Conventional wisdom is that increasing temperature causes quantum coherence to decrease. Using finite tem- 
perature perturbation theory and exact calculations for the strongly correlated bosonic Mott insulating state we 
show a practical counter-example that can be explored in optical lattice experiments: the short-range coherence 
of the Mott insulating phase can increase substantially with increasing temperature. We demonstrate that this 
phenomenon originates from thermally produced defects that can tunnel with ease. Since the near zero temper- 
ature coherence properties have been measured with high precision we expect these results to be verifiable in 
current experiments. 



PACS numbers: 05.30.Jp, 03.75.Hh, 03.75.Lm 



I. INTRODUCTION 

A key signature of the superfiuid to Mott insulator phase 
transition was the disappearance of long-range phase coher- 
ence, identified through the lossjrf sharp peaks in the inter- 
ference pattern after expansion JH]. Subsequent experiments 
showed HI, that a weak and slowly modulated inference 
pattern remained in the insulating regime, arising from short- 
range coherence between neighboring sites established by vir- 
tual hopping events [see Fig. Ufa)]. These processes occur 
with amplitude J/U, where J and U are the tunneling matrix 
element and onsite interaction energy, respectively [4]. Since 
J/U is small in the insulating regime, this tunneling process is 
well described by first order perturbation theory. Experiments 
in two-dimensional (2D) {Hi and three-dimensional (3D) J2l 
lattices have accurately measured this coherence and obtained 
quantitative agreement with the predicted J/U scaling. 

Another topic of much recent interest in optical lattice re- 
search has been the influence of temperature on the properties 
and phases of the many body system ]5|], with considerable 
development of techniques for measuring thermal properties 
JfJ. An area of particular interest has been the use of lattice 
loading |0] or other methods ^ to produce states of low tem- 
perature/entropy, e.g. to produce a system suitable for simu- 
lating spin Hamiltonians J^] . 

Here we study the influence of temperature on short-range 
coherence in the Mott insulating state. The physics we explore 
is a competition between two processes: (i) Gapped tunnel- 
ing, whereby an energy cost of U is needed for a particle to 
tunnel to its neighbor [Fig. Ufa)]. As mentioned above, this 
occurs with amplitude J/U and is the dominant process at 
T = 0. (ii) Gapless tunneling, occurs when the unperturbed 
state has a hole [Fig.QJb)] or particle [Fig.QJc)] defect which 
can tunnel with no energy cost and contribute strongly to co- 
herence. However, the production of such defects are expo- 
nentially suppressed at low temperatures. Here we study the 
Mott insulator coherence using finite temperature perturbation 
theory and exact solutions for small systems. We find that 
at T ~ O.lU/ks, gapless tunneling becomes important and 
can lead to an appreciable increase in the short-range coher- 
ence over the zero temperature value. We present results for 
both homogeneous (translationally invariant) and inhomoge- 
nous (harmonically trapped) lattices like in experiments y|]. 
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FIG. 1. (Color online) Schematic of tunneling processes that lead 
to short-range coherence being established, (a) Gapped tunneling: 
Particle-hole formation caused by tunneling, J, but costing the gap 
energy U. (b,c) Gapless tunneling: Thermally activated defects can 
tunnel without energy cost. 



H. THEORY 

We first consider a system of bosons in a transitionally in 
variant lattice described by the Bose-Hubbard Hamiltonian 

U 
~2 



H=-J^ 



a] a 



(i) 



We assume for simplicity that the lattice is simple cubic in ge- 
ometry and note that the tunneling is limited to nearest neigh- 
bor lattice sites. Our regime of interest is the Mott insulating 
phase (J/U <C 1) with an integer number (n) of atoms per 
site. To describe the short-range coherence between nearest 
neighbor sites, we evaluate the correlation function 



ci = (a] hi' 



(2) 



where i' is taken to be a nearest neighbor lattice site of i. 



A. Zero temperature 

We first review the zero temperature limit (see I2, [H [HI 
[TUl ). In this case, to lowest order in J, the manybody state 



2 



is a perfect Mott insulator with precisely n atoms per site, 
i.e. |V> (0) ) = IL [(a!)"/^] |vac), for which 



,(o) 



0. 



First order corrections to the ground state arise from the tun- 
neling of a particle between neighboring sites. This process 
is gapped (costing energy U to create a particle-hole excita- 
tion) and so the amplitude for the particle-hole admixture [see 
Fig-E a )] is \Jn(n + T)J/U (according first order perturba- 
tion theory), giving = 2n(h + 1) J/U . 

The effect of this neighboring site coherence is directly ob- 
servable in the momentum distribution as 



p(k) = N\w(t 



1 H cos(fc a a) 

Co , 

a— x,y,z 



(3) 



where N is the total number of atoms, cq = (a] Hi) (which is 
h in the Mott insulating regime), tu() is the Fourier transform 
of the Wannier state (the localized state at each lattice site), 
and a is the lattice constant. By fitting the momentum distri- 
bution ci was extracted in Ref. JH verifying the zero temper- 
ature prediction. Equivalently c\ can be determined from the 
visibility of the interference fringes (e.g. see ll2l [lol[l2ll ). 

Higher-order processes can establish coherence to next- 
nearest neighbors, leading to more rapid modulations of the 
momentum distribution [e.g. cos(2fc Q a)]. However, these 
contributions are much smaller and thus more difficult to mea- 
sure except for larger values of J/U where the system is near 
the superfiuid transition. For this reason we focus upon the 
nearest neighbor coherence c\. 



B. Finite temperature 
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FIG. 2. (Color online) Nearest neighbor coherence, (a) Perturbative 
treatment: (solid) full result , (dashed) gapped c{ a \ and (dash- 
dot) gapless contributions for the case of fj, = 1.5(7. (Black) 
show full perturbative calculation [Eqs. Q, {8}] while (yellow/grey) 
show the simple analytic results [Eqs. {9}, J 10H . (b) Variation with 
chemical potential and comparison of (solid) full perturbative results 
against (dots) exact calculations for a ID lattice, /i = 1.25(7 (black), 
1.5(7 (red/grey), and 1.75(7 (green/light grey). Inset: comparison of 
(solid) full perturbative results against (dots) exact calculations near 
the phase boundary /i = 1.95(7 (blue/gray), 1.05(7 (red/light grey). 
Other parameters: J = 0.01(7 and n — 2. 



At finite temperature the first order correlation function is 

* = (a]a t ,) = Z- x Tt{a\SLve-^->^}, (4) 

where Z = Tr{e^ l3i - H ^ tJ ' N ^} is the partition function, with 
(3 = 1/ksT, N = ^jOjCLj, and \i the chemical potential, 
used to adjust the average number of particles. 
Neglecting tunneling: By neglecting tunneling we obtain a 
system which can be solved exactly, but has no short-range 
coherence c\ = 0. However, this system is useful to consider 
as it is the unperturbed system for a perturbative expansion in 
the tunneling parameter. 

For J = the lattice sites decouple, and the problem re- 
duces to a single site anharmonic oscillator with energy levels 

E n = -jUn{n - 1) - fin, (5) 

for which the partition function is Z^ ' = ^2 n e~ l3Erl . The 
mean number of atoms per site is given by c = J2 n n Pn> 
where p„ = e~^ En /Z^ is the Gibbs factor. 
First order perturbative treatment: We use imaginary time 
perturbation theory to obtain an expansion for ©. At first 
order in J we have 

4 1 ' = jJ2n'(n + l)p nPn , f dTe-«+<~>\ (6) 



where n (n') is the number of atoms at site i (i'), and we 
have introduced the particle and hole energies given by e+ = 
E n+ i — E n , and e~ = E n -i — E n , respectively. Note that 
e n + e n> i s m e ener gy required (at zeroth order in J) to tunnel 
a particle from site i' to i, i.e. effect the change (n, n') — > 
(n + l,n' — 1). Eg. © is derived in a similar manner to 
expression (18) in lfl2ll and (44) in lfl3ll . and we do not repeat 
the details here. 

To evaluate the integral in Eq. (O we note the energy cost 
term is given by e+ + e~, = U(l + n — ??'), and thus two 
cases emerge: (a) Gapped excitations with n ^ n' — 1, and 
(b) gapless excitations for the case n = n' — 1. The first order 
coherence arising from the gapped excitations is given by 

(io) 2J ^ n'{n + l)p n p n > 

c i =77 l^—rr —(l-S n ,n'-i), (7) 

(7 * — ' 1 + n — n 

nn' 

while for the gapless excitations 

c[ lb) =(3 J £(n + lfp nPn+1 . (8) 

n 

In Fig. [2J a ) we show c\ and the separate contributions of the 
gapped and gapless parts, while in Fig. |2jb) we examine the 
dependence of c\ on chemical potential. The key result of the 
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paper, which we further analyze below, is that the gapless con- 
tribution at finite temperature serves to enhance system coher- 
ence over the T — state. For the various values of chemical 
potential we find that the maximum temperature induced co- 
herence occurs for T rj [0.1 — 0.2]U/kB, while the melting 
temperature for a trapped Mott insulator is ~ 0.2U /ks Q- 

We note that, by not explicitly accounting for the gapless 
excitations, the first order coherence expressions developed 
in lfl2tl and lfl3ll are formally divergent in the finite tempera- 
ture regime. Indeed, results for the finite temperature visibility 
presented in Fig. 4 of lfl2tl do not show the enhancement we 
observe at finite temperature, and instead show the character- 
istic behavior of the gapless part (c^ la - ) ) only. 

Since the first order expression for the coherence involves 
only a single tunneling event, it is independent of lattice di- 
mension. Higher order processes can explore the various pos- 
sible paths to tunnel between lattice sites, which is strongly 
dependent on the lattice dimension and geometry. However, 
the next correction to c\ is at 3rd order and for low tempera- 
tures and small J/U this should be negligible. 
Simple analytic treatment: Here we present an analytic sim- 
plification of results (0 and (© that helps clarify the under- 
lying physics. In the low temperature limit, ksT <C U, the 
summation in is dominated by the n = n' = n term 
(since p ¥l as 1) and we have 



2J , ^2 
— n{n + l)p n 



(9) 



which reduces to the known zero temperature result as T —> 0. 
The summation in c9~ b ^ is dominated by the n = n' ± 1 = n 
terms, for which pn±i ~ e _/3efi and we have 



Jib) 



f3Jpn{(n 



I) 2 



Pn+l 



n 2 pn- 



(10) 



Due to the energetic degeneracy of the gapless states coupled, 
this term occurs with a prefactor of /3 J rather than J/U, how- 
ever because defect states (which allow gapless tunneling) are 
thermally generated, the are exponentially suppressed by the 
factor Pn±i with temperature, and hence cj lb ' — > as T — >• 0. 

In Fig. 0a) we compare the simplified analytic results to the 
perturbative calculations. To do this we have taken p n = 1/z, 
p. n ±i = e^ 13 ^™ / z, where z = 1 + e _/3e ^ + e - ' 3 ' 5 ™ . Agreement 
is quantitatively good up until T ~ Q.2U /ks, as above this 
temperature the pn±i terms are becoming comparable to pn, 
and other terms neglected in our simplified treatment become 
important as the system begins to melt. 

The n atom-per-site Mott state occurs for chemical poten- 
tials in the range p/U <G (fi—1, n). The results in Fig. |3J a ) are 
evaluated in the middle of the n = 2 range. It is of interest to 
understand how the coherence varies with chemical potential 
for application to the trapped case (see Sec. Ill Cl . In Fig. [2b) 
we consider the behavior of cj 1 for several values of chemical 
potential. These results show that the chemical potential can 
considerably enhance the finite temperature coherence peak. 
This is easily understood as changing the chemical potential 
makes it easier to produce defects by altering the values of 
e^, e.g., for the case of p = 1.751/ we have that = 0.25C7 



(c.f. = 0.75U), leading to an enhancement of the first 
term in Eq. ( [Tol l. We note that the perturbative expansion is 
not valid when p/U is nearly integer, and in the trapped sys- 
tem so-called superfluid shells form where the local chemical 
potential satisfies this condition. 

Nonperturbative treatment: We also solve Eq. © by ex- 
act diagonalization for small ID systems in a lattice with 
M sites and periodic boundary conditions. In practice we 
do this using the truncated Fock basis |m, ...,tim) with 
< rij < n max to represent the grand canonical Hamilto- 
nian which we then numerically diagonalize. We retain the 
Ne lowest energy eigenvalues and vectors, {£ a , |^ a )}, and 
construct the density matrix p = Z^ 1 Yl a Zi e~^ £a \ip a ) (i/j a \, 
with Z = J2a=i e~ l3£a , from which we calculate c\. There 
are two important criteria for the validity of this calculation to 
describe c\ in the uniform (infinite) ID system: (i) M ^> 1, 
so that finite size effects are negligible for the calculation of 
ci O; (ii) The basis is sufficiently large to include all en- 
ergy scales accessed at finite temperature, i.e. [U n max (n max — 
\)/2,E Ne -£i] » k B T. 

In Fig.|2b) we compare the perturbation theory against the 
analytic results using AI = 5, n max = 7, and Ne = 620 for 
which £n e — £\ > 5U, We find good agreement between 
the exact and perturbation results up to the maximum tem- 
peratures considered, with the differences always remaining 
less than 3%, demonstrating that for these values of tunnel- 
ing the system is well within the perturbative regime. We find 
that perturbation theory starts to breakdown for J/U > 0.1. 
The inset of Fig. EJb) considers cases where p/U approaches 
an integer value. This shows that perturbation theory breaks 
down near the phase boundary for low temperature, and care 
must be taken to exclude such regions in the extension of the 
perturbation theory to a inhomogeneous system ifloll . While 
our comparisons here are for a one-dimensional (ID) lattice, 
it is worth noting that the agreement between perturbative and 
exact treatments is expected to improve with increasing lattice 
dimension IU3I1. 



C. Inhomogeneous system 

In experiment an additional external potential V(r) is 
present due to magnetic confinement and/or the Gaussian pro- 
file of the lattice beams. This can be represented in our the- 
ory by the replacement pN — > J^j Mi&jA; m ®> wner e 
pj = p — V(rj) is the local chemical potential, with V(rj) 
the value of the external potential at lattice site j. The gener- 
alization of Eq. (0 is then 



p(k) = N\wQf 



E 



2Ci„ 
N 



■ cos(fc Q a) 



(11) 



where C\ ja = Y^/w) W^i')> w i tn indicating a sum 

over nearest neighbor lattice sites along the a direction. In 
experiments lj]-[3|] the external potential varies slowly across 
the system, which allows us to use the local density approxi- 
mation. Thus we have Ci, Q = 2 f_ dp r g(p r )ci(p r ), where 



4 




across fi r ]. From Fig. [3] we can see that the enhancement of 
short-range coherence is around ~ 15% of the T = value 
for the 2D case and ~ 28% in the ID system. This could be 
easily observable with current experiments, particularly with 
the recent experimental improvements reported in 1 1611 . 



III. CONCLUSIONS 
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FIG. 3. (Color online) Nearest neighbor coherence in the x direc- 
tion. (Black dots) Homogeneous case [independent of dimension] 
with fi = 0.8(7. Inhomogeneous ID case for /i = 0.8(7 (green/light 
gray line), and for /i = 0.5(7 (blue/black dotted line). Inhomoge- 
neous 2D case for n = 0.8(7 (red/grey dashed line). Other param- 
eters: J/U = 8.35 x 10 -3 , fi values chosen so that n — 1 at trap 
center, interaction parameters are for 87 Rb in a lattice produced using 
counter-propagating A = 820 nm lasers, and harmonic confinement 
frequency is u = 2n x 40 Hz. 



g(Hr) is density of trap-potential states. We present results for 
C[ 1] x /N (calculated using first order perturbation theory lfl5ll ) 
in various systems in Fig. [3] where we notice a decrease in the 
peak value for the inhomogeneous case. This can be explained 
by the tendency of the spatially varying chemical potential to 
smear out the temperature induced coherence [e.g. consider 
the curves for various /i values in Fig.|2b)]. This smearing 
is least pronounced for a ID system as the density of states 
is strongly peaked around the chemical potential at the trap 
center [in contrast to 2D where the g(n r ) is evenly distributed 



In this paper we have studied the role of temperature on the 
coherence of a strongly correlated quantum system. Our main 
tool, finite temperature perturbation theory, provides insight 
into the underlying physical origin of the increase in coher- 
ence at finite temperature: gapless tunneling of thermally ac- 
tivated defects. We have verified our predictions against exact 
numerical calculations. While the formalism for finite temper- 
ature perturbation theory has been developed by several other 
authors, it is usually applied at T = with gapped excitations. 
The only application to finite temperature regime (that we are 
aware of) was a calculation for c\ in Ref. 11211 . Their results 
disagree with our work (predicting c\ to decrease with T) due 
to neglect of gapless excitations. 

Our study shows that the equilibrium state coherence of 
a Mott insulator state increases with temperature, propor- 
tional to the number of thermally induced defects. Interest- 
ingly, if defects are produced in the system by other means 
e.g. through non-adiabatic changes of the lattice potential, or 
three-body loss, then additional enhancement of the coherence 
is expected. Thus coherence measurements might comple- 
ment in situ imaging II 1711 to characterize defects in experi- 
ments. 
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